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Abstract 

In this paper, we calculated the B — ► (J/^f,rj c )K decays in the perturbative QCD approach 
with the inclusion of the partial next-to-leading order (NLO) contributions. We found that 
(a) when the large enhancements from the known NLO contributions are taken into account, 
the NLO pQCD predictions for the branching ratios are the following: Br(B° — > J/^K°) = 
5.2+H x 10~ 4 , Br(B+ -> J/^K + ) = 5.6 +3 2 7 g x 10~ 4 , Br(B° -» rj c K°) = 5.5±l;g x 10" 4 , 
Br(B + -> n c K + ) = 5.9±|;f x 10" 4 , which are roughly 40% smaller than the measured values, 
but basically agree with the data within 2 — a errors; (b) the NLO pQCD predictions for the 
CP-violating asymmetries of B — > (J/^,r] c )K decays agree perfectly with the data. 
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I. INTRODUCTION 



The B — > J/^K and B —>■ rj c K decays are phenomenologically very interesting decay 
modes and have drawn a great attention for many years. Although the underlaying weak 
decay of b — > ccs is simple, but a clear understanding of the exclusive B — ► (J/^,1] C )K 
decays is really difficult because of the involving of the complex strong-interaction effects. 

On the experiments side, the experimental studies for the "Golden-plated" B — > 
J/^Ks^l decays result in the precision measurement of sin2/?[l[. The branching ratios 
of B —>■ (J/ty,rj c )K decays and other similar decays involving a charmonium and a light 
pseudo-scalar or vector meson as the two final state meson, have been measured with 
good or high precision [j], 0]: 

Br(B° -> J/^K°) = (8.71 ± 0.32) x 10" 4 , 

Br(B + -> J/^K + ) = (10.07 ±0.35) x 10" 4 , (1) 

Br(B° -> r] c K°) = (8.9 ± 1.6) x 10" 4 , 

Br(B + -> 7^+) = (9.1 ± 1.3) x 10" 4 . (2) 

The accuracy of above measurements will be improved rapidly along with the running of 
the relevant LHC experiments. 

On the theory side, such B meson charmonia decays have been studied intensively by 
employing various theoretical methods or approaches, for example, in Refs. [!, 0, 0, @, 
0, IE IE 12; EH 12, 13l . But unfortunately, it is still very difficult to give an satisfactory 



explanation for the corresponding data without the worry about the serious problems. 

For B — > J/^/K decay, for example, the theoretical predictions for its branching ratio 
in both the naive factorization approach (NFA) jij and the QCD factorization (QCDF) 
approach (l4| are much smaller (a factor of 7 ~ 10 ) than the measured values^: BriB — > 
J/^K) ~ 1.1 x 10~ 4 when the twist- 2 distribution function (DA) 4>j/^(x) = 6x(l — x) 
was employed [5j. 

In Ref. [6|, the authors studied the effects of twist-3 DA cf>^ and found that the resultant 
enhancement to the Wilson coefficient d2{J/^K) and consequently to the branching ratio 
BriB — > J/^/K), induced through the spectator diagram, can be large. But one should 
note that there are also logarithmic divergences arising from spectator interactions due 
to kaon twist-3 effects, this is always a serious problem in the QCDF approach. 

In Refs. 0, @|, the authors studied the decays B — > (j] c , rf c , Xco> Xci)K m the QCDF 
approach and found that (a) the logarithmic divergences will arise from the spectator 
interactions due to the kaon twist-3 effects; and (b) the predicted decay rate is Br(B — > 
rj c K) = 1.9 x 10~ 4 , which is still a factor of 5 smaller than the measured value in Eq.(T5]). 
They concluded that the QCDF approach with its present version can not be safely applied 
to exclusive decays of B meson into charmonia [9| . 

The B — > J/^/K decays have also been investigated by employing the QCD light- 
cone sum rules (LCSR) [lOj. The authors calculated the nonfactorizable contributions 
to the B — > J/^/K decay coming from the exchanges of the soft gluons between the 
emitted J/\l/ and the kaon. But their predictions for branching ratios is still too small, 
Br(B — > J/tyK) y 3.3 x 10~ 4 , to accommodate the data. 

In Refs. [I2I. 13| . the authors studied B — > ( J/\l/, r] c , Xco,ci)K^ decays in a formalism 



that combines the QCDF factorization and the perturbtive QCD (pQCD) approaches [111]. 
They employed the QCDF approach to calculate the factorizable contribution, but the 
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pQCD approach to evaluate the nonfactorizable corrections to the considered decays. Ac- 
cording to their studies 12, 13[ we see that (a) the theoretical predictions for the branching 



ratios of B — > (J/^/,Xco,ci)K decays can be large and consistent with the data; (b) the 
B — > r) c K decays still exhibit a puzzle: the predicted result is BriB — > r] c K) ~ 2.3 x 10 -4 , 
much smaller than the measured values as given in Eq. (pQ). Furthermore, it should be 



mentioned that these results [12j, [13j were obtained by treating one decay with two dif- 
ferent factorization approaches: the self-consistency of such "mixing-approach" may be a 
serious problem. 

Up to now, a clear and satisfactory theoretical interpretation for the measured large de- 
cay rates of B — > (J/^,rj c )K are still absent. We call this situation the U B — > (J/^, T] C )K" 
puzzle. In this paper, we will calculate the branching ratios and CP-violating asymme- 
tries of the four B — ► (J/^f, Vc)K decays by employing the pQCD factorization approach: 
(a) we evaluate both the factorizable and nonfactorizable contributions in the pQCD ap- 
proach; (b) besides the full leading order (LO) contributions in the pQCD approach, the 



currently known next-to-leading order (NLO) contributions 15| (specifically the QCD 
vertex corrections for the considered decays) are also included. 

The paper is organized as follows: in Sec. [TTl we firstly present the formalism of the 
pQCD approach, and then make the analytic calculations and show the decay amplitudes 
for the considered decays. In Sec. IIII[ we show the numerical results and compare them 
with the measured values. A short summary and some conclusions are given in the last 
section. 

II. FORMALISM AND PERTURB ATIVE CALCULATIONS 
A. Formalism 

In recent years, the pQCD factorization approach has been used frequently to calculate 
various B meson decay channels. For the two body charmless hadronic B meson decays 
the pQCD predictions for the branching ratios and CP- violating asymmetries generally 



agree well with the measured values [16|, ll7|, [18|, [19|, |20|, |2l|. In Ref. (22j, the authors 
calculated B -> D*K, Di* )+ Di* y and B s -> £>(*)+£>(*)- decays and found that the pQCD 



approach works well for such decays. In a previous paper [23j, the B — > J/^K decays 
have been studied by employing the pQCD approach at leading order. Here we try to 
apply the pQCD approach to calculate the B — > (J/^,rj c )K decays with the inclusion of 
the NLO corrections. 

In pQCD approach, the decay amplitude of B — > M2M3 decays 1 can be written con- 
ceptually as the convolution, 

A{B - M 2 M 3 ) ~ Jd A k l( l 4 k 2 d% Tr [C(t)$ B (fc 1 )$ M2 (A; 2 )$ A/3 (fc3)^(fci,fc 2 , h,t)} , (3) 

where the term "Tr" denotes the trace over Dirac and color indices. C(t) is the Wilson 
coefficient which results from the radiative corrections at short distance. In the above 



Here M 2 — (J/^, T]c) is the emitted charmonium, and M 3 is the kaon which absorbed the spectator 
quark. 
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convolution, C(t) includes the harder dynamics at larger scale than mg scale and describes 
the evolution of local 4-Fermi operators from myy (the W boson mass) down to t ~ 

O^Athb) scale, where A = tub — m&. The function H(ki, k 2 , k 3 , t) is the hard part and 
can be calculated perturbatively. The function $m is the wave function which describes 
hadronization of the quark and anti-quark to the meson M. While the function H depends 
on the process considered, the wave function $m is independent of the specific process. 
Using the wave functions determined from other well measured processes, one can make 
quantitative predictions here. 

Using the light-cone coordinates the B meson and the two final state meson momenta 
can be written as 

P 1 = ^|(l,l,O r ), P 2 = ^|(l,r 2 ,0 T ), P 3 = ™|(0,l-r 2 ,O T ), (4) 

respectively, where r = mM 2 / m B, and the light pseudo-scalar meson masses mM 3 = mx 
have been neglected. The longitudinal polarization of vector J/^f, 6l, is given by = 
— (1, — r 2 ,^, Ot). Putting the light (anti-) quark momenta in B and M 3 mesons as 

k\ and k 3 , respectively, we can choose 

h = (x!P+, 0, k ir ), k 3 = (0, x 3 P 3 , kgr). (5) 

For M2, the momentum fraction of c quark is chosen as X2P2. Then, the integration over 
ki, &2, and k 3 in Eq.Q will lead to 



A(B — > M 2 M 3 ) ~ J dxidx2dx 3 bidbib 2 db2b 3 db 3 

•Tr [C(t)$ B ( Xl , 6i)$ A / 2 (x 2 , b 2 )^M 3 (x 3 , b 3 )H(x h b u t)S t ( Xi ) e' 8 ^] ,(6) 

where 6, is the conjugate space coordinate of k iT , and t is the largest energy scale in 
function H(xi, fej, t). The large logarithms \n(mw/t) are included in the Wilson coefficients 
C(t). The large double logarithms (ln 2 Xj) on the longitudinal direction are summed by 
the threshold resummation [24j, and they lead to S t (xi) which smears the end-point 
singularities on Xj. The last term, e~ s ^\ is the Sudakov form factor which suppresses the 
soft dynamics effectively [25|. Thus it makes the perturbative calculation of the hard 



part H applicable at intermediate scale, i.e., m B scale. 

For the considered decays, the weak effective Hamiltonian Ti e ff for b — > s transition 
can be written as 



n, GF 

Heff ~ 7! 



10 



v: b V cs (CMOt(fi) + Cfa)Ofoi)) - V*V ts Ciin) Oi(ii) 



i=3 



(7) 



where Ci(n) are Wilson coefficients at the renormalization scale /1 and Oj are the four- 
fermion operators: 



0{ = 


Scfl^Lcp ■ 


cpl^Lba , 

Eg' q'plM'p > 


C 2 


= s a ^Lc a ■ cprfpHjp , 


o 3 = 


s a ^Lb a ■ 


4 


= SoTf^Lbp ■ E g ' q'piM* > 

= sa^Lbp ■ E g ' q'p^Rq' a , 


5 = 


Sa^Lba ■ 


o e 


o 7 = 


ls a ^Lb a 


■ Eg' eq'Q'piMp , 


o 8 


= \s a ^Lb p ■ E g / e q ,q'^Rq' a 


o 9 = 


\s a ^Lb a 


O10 
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FIG. 1: Typical Feynman diagrams contributing to B — * (J/^ ,r/ c )K decays at leading order. 



where a and (3 are the SU(3) color indices; L and R are the left- and right-handed 
projection operators with L = (1 — 75), R = (1 + 75). The sum over q' runs over the quark 
fields that are active at the scale \i = 0(m b ), i.e., q'e{u, d, s, c, b}. 

In PQCD approach, the scale "t" appeared in the Wilson coefficients Cj(t), the hard- 
kernel H(xi, bi, t) and the Sudakov factor e~ s ^ is chosen as the largest energy scale in the 
gluon and/or the quark propagators of a given Feynman diagram, in order to suppress 
the higher order corrections and improve the reliability of the perturbative calculation. 
Here, the scale "i" may be larger or smaller than the scale. In the range of t < nib or 
t > mb, the number of active quarks is Nf = 4 or Nf = 5, respectively. For the Wilson 
coefficients Cjjix) and their renormalization group (RG) running, they are known at NLO 
level currently [26J. The explicit expressions of the LO and NLO Cj(ms) can be found 
easily, for example, in Refs. j26|, 13 . 

When the pQCD approach at leading-order are employed, the leading order Wilson 
coefficients Ci(m\y), the leading order RG evolution matrix U(t,m)( ' from the high scale 
m down to t < m and the leading order a s (t) are used: 



a s (t) 



An 



A) In [tV^Qco] 



(9) 



where (3 = (33 - 2JV»/3. 

When the NLO contributions are taken into account, however, the NLO Wilson coef- 
ficients Ci{mw)i the NLO RG evolution matrix U(t,m,a) ( see Eq. (7.22) in Ref. [2o| ) 
and the a s (t) at two-loop level will be used: 



a s (t) 



An 



A,ln[t 2 /^cd 



Pi l n[ln[t 2 /^ cc ]] 



where /3 = (33 - 2JV»/3, A = (306 - 38JV»/3. By assuming A^ CD 
will get K { q> cd = 0.287 GeV (0.326 GeV) for LO (NLO) case. 



(10) 



0.225 GeV, we 



As discussed in Ref. 20], it is reasonable to choose fio = 1-0 GeV as the lower cut-off 
of the hard scale t. In the numerical integrations we will fix the values Cj(t) at Cj(1.0) 



whenever the scale t runs below the scale /x = 1-0 GeV [2CJ, [21|, unless otherwise stated. 



B. B —> J/^Mz decays at leading order 

At the leading order pQCD approach, as illustrated in Fig. 1, the relevant Feynman 
diagrams for the considered decays include the factorizable emission diagrams (Figs, la 
and lb) and the non-factorizable spectator ones (Figs.lc and Id). The operators 0\^, 
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03,4 and 09,10 are the (V — A)(V —A) currents, while Os,6 and O73 are the (V — A)(V+A) 
currents. By analytic calculations of Fig. la and lb, one finds the corresponding decay 
amplitudes 

/l ^oo 
dx x dx 3 J M&iM& 3 0b(^iA) 

[(1 - r 2 )(l + x 3 ) - x 3 r 2 ] <\> A Mz [x 3 ) + r (l - 2x 3 ) [<p P Mz {x 3 ) + <&(x 3 )] 



-r r 2 [(1 - 2x 3 )<fa a (x 3 ) - (1 + 2x 3 )<j ) T M3 (x 3 )} 
■a s {t\) h e (x 1 ,x 3 , bx, b 3 ) exp [S ab (tl)] 
+2r (l - r 2 ) 0m 3 (x 3 ) • a s (tg)/i e (x 3 ,xi,6 3 ,6i)exp [S ab (t 2 e )] L (11) 



where r = /rug, and 0^ = 4/3 is a color factor. The hard function h e , the scales t l e 
and the Sudakov factors S a b are displayed in Appendix lAl 

Now we consider the contributions of the operators 05,6,7,8 i n the Fig.l. In some decay 
channels, some of these operators contribute to the decay amplitude in a factorizable way. 
Since only the vector part of (V + A) current contribute to the vector meson production, 
(M 3 \V- A\B)(J/V\V + A\0) = (M 3 \V- A\B)(J/V\V - A\0), that is 

f v+a _ f v-a / 12 n 

For the non-factorizable diagrams 1(c) and 1(d), all three meson wave functions are 
involved. The integration of b 3 can be performed using 6 function 5(b 3 — £>i), leaving only 
integration of b\ and b 2 . For the (V — A)(V — A) operators, the corresponding decay 
amplitude is 



16y/E 



M j/<km 3 = ^~ nC F m B / dx 1 dx 2 dx 3 \ b 1 dbib 2 db 2 (j) B (x 1 ,b 1/ 



jo 

A („ \ a„„ „ it („ \ jT 



x|2rr c J/ ^(x 2 )0 M3 (x 3 ) - 4rr o r c J/v j,(x 2 ) 4>m 3 { x z) 

- [x 3 + 2(x 2 - x 3 )r 2 ] (j)j/q,(x 2 )(f)M s (x 3 ) 

+2r [x 3 + (2x 2 - x 3 )r 2 ] (t) L j/^{x 2 )(t) T M3 (x 3 ) 

■a s (tf)h f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 )exp[-S cd (t f )] , (13) 

where r c = m c /niB and m c is the mass for c quark. 

For some decay channels, the (S — P)(S + P) operators can be obtained from the 
(V — A)(V + A) operators by making the Fierz transformation, in order to get right color 
and flavor structure for factorization to work. For these (S — P)(S + P) operators, the 
corresponding decay amplitude can be written as 

M s+P - -M V ~ A (14) 
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For B — > J/^>M 3 decays, by combining the contributions from different Feynman 
diagrams, the total decay amplitude can be written as 



M(B - J/*M 3 ) = Fy- A M Jj/y {V: b V cs a 2 - V*V ts (a 3 + a 5 + a 7 + a 9 )} 



+m; 7 -^ 3 {v; b v cs c 2 - v; b v ts (c 4 - c 6 - c 8 + c 10 )} , (is) 

where a,- is the combination of the Wilson coefficients C,: 



C2 C 1 

a 2 = Ci + — ; o» = Ci + — jp» /or i = 3,5,7,9, 

where C 2 ~ 1 is the largest one among all Wilson coefficients. 



(16) 



C. B — > ry c M3 decays at leading order 

Following the same procedure as for — > J/^/K decays, it is straightforward to cal- 
culate the decay amplitudes for B — > r/ c M 3 decays. 

M(B - 77 C M 3 ) = F^/, c [V^ a 2 - V^V ta (a 3 + a 5 + a 7 + a 9 )] 

+<Ma [KtK s c 2 - v; b v ts (c 4 + c 6 + c 8 + do)] , (17) 

where the functions F^j£, M^J^, etc, are of the form 



f v-a 

ricM 3 



"1 

—FYt^ = 87rCFm 2 B / dx\dx 3 I bidbib 3 db 3 (f) B (xi,bi) 



ricM 3 



X 



[(1 _ r 2) (1 + X3 ) _ X3r 2] ^ (xs) + r()(1 _ 2x3) [0^^) + 0^(x 3 )] 



+r r 2 [(1 + 2x 3 )0^ 3 (x 3 ) - (1 - 2x 3 )0^ 3 (x 3 )] 
•a s (tg) h e (zi, ar 3 , b u b 3 ) exp[-S ab (tl)] 
+2r (l -r 2 ) 0^3(^3) • a s (t 2 )/i e (x 3 ,a;i,6 3 ,6i)exp[-S' a6 (t 2 )] L 



(18) 



m v-a _ m s+p _ _ We 



pi pco 

irC F m B / dxidx 2 dx 3 / b 1 dbib 2 db 2 <pB{xi,bi)4> v rjc (x 2 ,b 2 ) 
Jo Jo 



xx 3 



(l-2r 2 ) 0^( X3 )-2r o (l-r 2 ) 0^(x 3 ) 



•oi s {tf)hf(x u x 2 , x 3 , 61, 6 2 ) exp[-5' cd (i/)] , 



(19) 



where is the leading twist- 2 part of the distribution amplitude for the pseudo-scalar 
meson i] c . 
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FIG. 2: Typical Feynman diagrams contributing to B — > M2M3 decays at NLO level. 



D. NLO contributions in pQCD approach 

For a gen eral B — > M2M3 decays, the power counting in the pQCD factorization 
approach [15J is different from that in the QCD factorization 14j. In the pQCD approach, 
the NLO contributions may include the following parts [15l. 



2. 



4. 



The Wilson coefficients C^m^) and the renormalization group evolution matrix 
U (£, m, a) at the NLO level, and the a s (t) at the two- loop level [26|] should be used. 

All the Feynman diagrams, which lead to the decay amplitudes proportional to 
should be considered. 

Currently known NLO contributions: (a) the vertex corrections; (b) the contribu- 
tions from the quark-loops and the chromo-magnetic penguins (0 8g ), as illustrated 
in Fig. F3J 

The NLO contributions can also come from the Feynman diagrams as shown in 
the Figs. 5-7 in Ref. |20J. The analytical calculations for these (more than 100!) 
Feynman diagrams have not been completed yet. 



For the considered B — > J/^/K and rj c K decays, only the vertex corrections (see Fig. 2a- 
2d) among the known NLO contributions will contribute. For the four vertex correction 
diagrams Fig.2a-2d, as was confirmed in Ref.Q], the infrared divergences from the soft 
gluons and collinear gluons in the four diagrams will be canceled each other, respectively. 
So the total contributions of these four figures are infrared finite. In other words, these 
vertex corrections can be calculated without considering the transverse momentum effects 
of the quark at the end-point region in collinear factorization theorem. Therefore, there is 
no need to employ the factorization theorem here. The vertex corrections to the B — > 
J/ipK decays, denoted as // in QCDF, have been calculated in the NDR scheme 0], 
and can be adopted directly. Their effects can be combined into the Wilson coefficients 
associated with the factorizable contributions: 



C 2 a s C F 

a2 = Cl + T r + ^N: C2 



48 + 12 \n^ + fj 



(20) 
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[1 (GeV) 



FIG. 3: Dependence of a 2 on the renormalization scale /t for (a) B — > J/^K and (b) B — ► ?7 C .K" 
decays. The solid ( dashed ) curve stands for a 2 at NLO (LO) level without the vertex corrections, 
while the dotted ( dash-dotted ) curve refers to the real (imaginary) part of the a 2 at NLO level 
with the vertex corrections. 



°- c =+| + sf c ^- 121 V-'')' (21) 

C$ a s C F f , mi, f 

a 9 = C 9 + ^ + ^C 10 f-18 + 121n^ + M , (22) 
N c 47r N c \ /j, 



with the function f J 



2V2K r t 

J 1 = j J dx 2 (j)j / y{x 2 ) 



3(1-2x2) , 2 , , 2rl(l-x 2 ) 

In x 2 — 3m + 3 ln(l — r 2 H ^ 

1 — x 2 1 — r^x 2 



■ (23) 



where r2 = mj/^/mB and those terms proportional to have been neglected. In Eqs.fl20I 
E2l . the Wilson coefficients Cj at NLO level should be used when the NLO vertex correc- 
tions are taken into account. 

For B — > f] c K decays, it is easy to obtain the corresponding NLO vertex corrections 
from those Wilson coefficients in Eqs. (I20H22I) . by the replacement of the parameter fj/^, 
and r 2 {J /^) in // with f Vc and r 2 (rj c ) |9j, respectively. 

Since both B — > J/^/K and r\ c K are color-suppressed decays, the Wilson coefficient 
a 2 in general plays the dominate role. It is instructive to check the variation of a 2 at 
LO or NLO level, with or without the inclusion of the NLO vertex corrections. From 
Figs. 3a and 3b, one can see that (a) the //-dependence of a 2 (J/^/K) and a 2 {j] c K) are very 
similar; (b) the //-dependence of a 2 is decreased effectively due to the inclusion of the 
NLO vertex corrections; and (c) the vertex correction provides a large imaginary part to 
a 2 , and therefore an effective enhancement of 1 0,2) due to the vertex correction is expected. 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



A. Branching ratios 

The input parameters and the wave functions to be used in the numerical calculations 
are given in Appendix [B] 

Firstly, we find the pQCD predictions for the corresponding form factors at zero mo- 
mentum transfer: 

F B - K (q 2 = 0) = 0.32^K) , (24) 

for fs = 0.19 GeV, and Ub = 0.40 ± 0.04GeV. It agrees very well with that obtained in 
QCD sum rule calculations j28| . 

Now we calculate the branching ratios for those considered decay modes. With the 
complete decay amplitudes, we can obtain the decay width for the considered decays, 

F(B -> M 2 K) = %^(1 - r 2 ) \M(B -> M 2 K)\ 2 , (25) 

where r = mj/^/rriB or m Vc /mB- 

By using the input parameters and wave functions as given in Appendix (Bj we find the 
LO and NLO pQCD predictions (in unit of 10 -4 ) for the CP-averaged branching ratios 
of the four B — ► J/^K and r] c K decays and show them in the Table [B The predictions 
listed in the column one are the LO pQCD predictions by setting fi = 1.0 GeV. The 
predictions listed in the column two is the NLO pQCD prediction by setting /i = 1.0. 
The first theoretical error in these entries arises from the B meson wave function shape 
parameter ut = 0.40 ± 0.04 and the decay constants fj/y = 0.405 ± 0.014 GeV and/or 
f Vc = 0.420 ±0.050 GeV. The second error is from the combination of the uncertainties of 
Gegenbauer moments af = 0.17 ±0.17 and/or a 2 = 0.115 ±0.115. In the fourth column 
of Table | as a comparison, we also cite the typical theoretical predictions obtained 
previously by using various approaches or models [1, @, 0, [l(| • 

In the last column of Table [TJ we list the world averages of the experimental mea- 
surements [H, 0|. One can see that the LO pQCD predictions for the branching ratios 
are indeed much smaller than the measured values, but the pQCD predictions with the 
inclusion of the vertex corrections can enhance the branching ratios evidently for both 
B — > J/^/K and B — > 7] C K decays. When the NLO enhancement are included, the pQCD 
predictions are basically consistent with the data within the still large theoretical errors. 
Of course, the central values are still about 40% smaller than the measured values. There 
is still some room left for non-perturbative contributions. 

Now we investigate in more detail why the vertex corrections can provide a significant 
enhancement. The total decay amplitudes as given in Eqs. f|T5|) and (fT7|) can be in re- 
written in the following form 

M = [F c + F P ] Fac + [M T + M P ] 5pec . , (26) 

where Fq and Fp stands for the "color-suppressed" and the penguin part of the factor- 
izable contribution, coming from the emission diagram Fig. la and lb; while My and Mp 
stands for the "Tree" and the penguin part of the nonfactorizable contribution, coming 
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TABLE I: The LO and NLO pQCD predictions (in unit of 1(T 4 ) of Br(B -» J/WK) and 
Br(B — > rj c K). For comparison, we also cite the typical theoretical predictions as given in 
previous literatures (in the fourth column), and the measured values [ll,0]. 



Channels 


LO 


NLO 


Others 


Data 


B° -> J/VK° 
B + -» J/VK+ 


i 1+0.8+0.9 
i - i -0.5-0.5 
i 9 +0.9+0.9 
i - z -0.5-0.5 


t- 9+1.0+3.4 
,J - z -0.9-2.6 
c c+1.0+3.6 
d -°-0.9-2.8 


~ 1.0 [5] 
~ 3.3 [10] 


8.71 ±0.32 
10.07 ±0.35 


B° -> 77 CJ ftT° 
B + -> r? c iT + 


n o+0. 5+0.1 
u -°-0.3-0.1 
n O+0.5+0.2 
u -°-0.2-0.1 


cr k+2.1+1.0 
°'°-1.7-1.0 
r q+2.2+1.2 


~ 2[9, L2] 
~ 2[10,i2j 


8.9 ± 1.6 
9.1 ± 1.3 



from the spectator diagram Fig.lc and Id. From Eq. (1151) . for example, it ie easy to 
separate the total decay amplitude A4(B — > J/^/K) into the following four parts 

Fc = FjfqrKfj/v V cb v cs a 2 ; 

F P = -Fj-^fj/* V t * b V ts (a 3 + a 5 + a 7 + a 9 ) ; (27) 

M T = M V j-„ A K V; b V cs C 2 ; 

M P = -Mj f ~ A K V t * b V ts (C 4 - C 6 - C7 8 + Cio) (28) 

By numerical calculations, we find easily the numerical values (in unit of 1CT 3 ) of the 
individual parts and the total decay amplitude of Ai(B — > J/^/K) at the LO and NLO 
level: 

M LO = -1.147 -0.092+ (1.487 + i0.529) + (0.046 + Z0.004), 



F C f p M T Mp 

= 0.294 + i0.533, (29) 
m nlo = (_ .546 - il.433) + (0.004 - i0.038) + (1.367 + i0.485) + (0.037 + z0.002), 

V «, ' V «, ' S «, ' S v ' 

F c Fp M T M P 

= 0.862 -^0.984, (30) 
and the ratio of the square of the decay amplitude J\A NLO and M. LO is 

\M NLO \ 2 

R M {B -> J/^K) = I— — j = 4.62. (31) 

For B — > r] c K decay, we find the similar result: Rm(B — > r) c K) = 7.0. 
From the above numerical results, it is easy to see that 

• As generally expected, both the the factorizable penguin contribution ( Fp) and the 
nonfactorizable part ( Mp) to the total decay amplitude are always small in mag- 
nitude ( less than 10%), when compared with the "Tree" and "Color-suppressed" 
parts (My and Fc)- 
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At the leading order, Fc = —1.147 is large in size, but largely canceled by the real 
part of Mt ( Re(Mr) = 1.487 ) when one sums up the factorizable and nonfactoriz- 
able contributions. This strong cancelation results in a small LO pQCD prediction 
for the decay rates. 

At the next-to-leading order, the NLO Wilson coefficients will be used. The penguin 
part Fp and Mp remain small in magnitude. The variations of Mt and Mp due 
to the replacement of the LO Wilson coefficients by the NLO ones are also small 
as expected 2 . For the "color-suppressed" Fc part, however, things become much 
different. The real part of Fc changes from —1.147 to —0.546, the previous large 
cancelation between the real parts of Fc and Mt become weak significantly, while a 
large imaginary part lm(Fc) = —1.433 is also produced. These two changes lead to 
a large lA'! 7 ^ ] 2 and consequently a large NLO pQCD prediction of the branching 
ratios. 

Although the NLO pQCD predictions for the branching ratios of the considered 
decays are consistent with the date within the still large theoretical uncertainties, 
but the central values of the NLO pQCD predictions, as listed in Table HJ are still 
about 60% of the measured values. Certainly, there are still some room left for the 
non-perturbative long distance effects or other unknown high order corrections. 

Among the three kinds of known NLO contributions in the pQCD approach, only the 
vertex corrections are relevant to B — > (J/^/,i] c )K decays and taken into account 
here. Other possible NLO contributions coming from the Feynman diagrams as 



shown in Figs. 5-7 in Ref. |20j are still unknown at present. But they are generally 
expected to be the small part of the NLO contributions in the pQCD factorization 
approach [TBI. [20| . 



In Ref. [29|, the authors attempted to estimate the soft rescattering effect for 
B + — > J/^>K + decay and concluded that such effect is comparable in size with 
the experimental one. But it is worth mentioning that the estimation as presented 
in Ref. [29] has a very large theoretical uncertainty. We believe that the long dis- 
tance effects for the considered decays are exist and may be large, but much more 
careful studies should be made before one can find an reliable estimation about 
them 0, Q E ■ 



As discussed in previous section, we choose /x = 1-0 GeV as the lower cut-off of the 
hard scale t. It should be noted that in the considered decay channels, the characteristic 
hard energy scale maybe smaller than 1 GeV, while all the meson distribution amplitudes 
are defined at 1 GeV. In the numerical integration, we therefore should set the Wilson 
coefficients Cj(i) = Cj(l GeV) whenever t < 1 GeV to ensure the reliability of our 
perturbative calculations. 

Now we turn to study the ratios of the branching ratios for some phenomenologically 
relevant decay modes. One advantage of the ratios is the possible cancelation of the theo- 
retical uncertainties of individual calculations. Using the data of the measured branching 



2 In the evaluation of M LO and M NLO , the LO and NLO Wilson coefficients C 2 and C 4)6 , 8 ,io will be 
used, respectively. 
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ratios as given in Ref. [H, the three ratios R^ xp can be defined as the following 

R T P = Tr-TTTT TTTTTk = 0-90 ± 0.13, 



Br(B + 


-> VcK + ) 


Br{B+ - 


* J/^K+) 


Br(B° - 


-> ^°) 


Br(B° - 


J/^K°) 


Br(B° - 


* VcK°) 



RT P = T~~7~~~77 J'l.l = 1.02 ±0.19, 



= ^1A±: ^£iL± = 0.88 ± 0.16 . (32) 

3 Br(B~ -> r] c K-) v 7 



TABLE II: The ratios Ri of the pQCD predictions for the branching ratios. The ratios as given 
in PDG 2008 [l| are also shown in the last column. 



Ratios 


LO 


NLO 


Data 


Ri 

R2 

Rs 


n 7 +0.9 

7+ 1 ' 3 

1 o +0 - 6 


l.UO_ 4 2 

i - UD -0 .43 
n OQ+0.52 
u.yo_ 27 


0.90 ±0.13 
1.02 ±0.19 
0.88 ±0.16 



By comparing the pQCD predictions of the ratios and the measured ones, as listed 
in Table [III we can see t na t (a) the consistency between the pQCD prediction for the 
ratios and the measured ones is improved significantly when the NLO contributions are 
taken into account; (b) the theoretical uncertainties are still large because of our poor 
knowledge about the Gegenbauer moments af 2 - 

In short, from the above pQCD predictions for the branching ratios and the detailed 
phenomenological analysis, we can conclude that the pQCD predictions for the branching 
ratios become close to the data due to the significant enhancement of the NLO vertex 
corrections. 



B. CP- violating asymmetries 



Now we turn to the evaluations of the CP- violating asymmetries of B — > M 2 M3 decays 
in pQCD approach. For the charged B meson decays, the direct CP- violating asymmetries 
A d Q P can be defined as usual. For both B + — > J/^>K + and t] c K + decays, there are no 
direct CP violation, since there is no weak phase appeared in their decay amplitude, as 
can be seen easily in Eqs. ( 1T5T) and (TT7]) . This theoretical expectation agrees well with the 
data HQ: 

Ac r p(B + -> J/^K + ) = 0.017 ± 0.016 , 

Ac r p(B+ -> r] c K + ) = -0.16 ± 0.08 . (33) 

which are consistent with zero within 2a errors. 

For the B° — ► M 2 M 3 decays, because these decays are neutral B meson decays, so we 
should consider the effects of B° — B° mixing. The direct and mixing induced CP- violating 
asymmetries Aq P and A^p can be written as 

A dir _ \^CP\ 2 ~ 1 .mix _ 2Im(X C p) , - 

Acp - i + \\ CP r cp ~i+\\ CP r { } 
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where the CP- violating parameter Xcp is 

_ VlV td (f\H eSS 0>) _ 2i9 (/|g e// |B°) 
^-'"^(/IffeZ/IB )-"' 6 ffUQS 5 )' <35) 

where ??/ is the CP-eigenvalue of the final states. By using the the input parameters as 
given in Appendix B, we find the following pQCD predictions 

A d S r P (B° - J/*K° S ) = A d S r P (B° - VcK° s ) « , 

-> J/^A"°) = -> 77c^) = (70.911?) %, (36) 



where the dominant error comes from p = 0.1351q;q^ and f\ = 0.349+q;q}7 pj. It is easy 
to see that the pQCD predictions for CP -violating asymmetries agree perfectly with the 
experimental measurements [H, @] ■ 



IV. SUMMARY 

In this paper, we calculated the branching ratios and CP-violating asymmetries of the 
four B — > (J/^/,r] c )K decays by employing the pQCD factorization approach with the 
inclusion of currently known NLO contributions. 

From our numerical calculations and phenomenological analysis, we found the following 
results: 

• The inclusion of the known NLO contributions can result in a factor of five enhance- 
ments to the leading order results. The NLO pQCD predictions for the branching 
ratios are the following 



Br(B° -» 


J/VK°) 


_ c 0+3.5 
— J.Z_ 28 


x 10" 


-4 


Br(B+ -> 


J/^K + ) 


— O.U_2.9 


x 10" 


-4 


Br(B° 


- VcK°) 


— 0.0_ 2 Q 


x 10" 


-4 


Br(B + 


- Vc K + ) 


= 5.91" 


x 10" 


-4 



Although the central values of the pQCD predictions are still 40% smaller than the 
measured ones, they basically agree with the data within 2a errors. One can also 
see that, on the other hand, there are still some room left for the non-perturbative 
contributions. 

• The pQCD predictions for the CP-violating asymmetries of the considered decays 
also agree perfectly with the data. 

• In this paper, only those currently known NLO contributions have been taken into 
account. To obtain a complete NLO calculations in the pQCD approach, the still 
missing pieces should be evaluated as soon as possible. 



Acknowledgments 

The authors are very grateful to Hsiang-nan Li, Cai-Dian Lii and Ying Li for valuable 
discussions. This work is partially supported by the National Natural Science Foundation 
of China under Grant No. 10575052, 10605012 and 10735080. 



14 



APPENDIX A: RELATED FUNCTIONS 



We show here the hard function hiS, coming from the Fourier transformations of the 
function H^°\ 

h e (x 1 ,x 3 ,b 1 ,b 3 ) = K (a/xix 3 (1 - r 2 )m B bi ) 9{bi - b 3 )K (a/x 3 (1 - r 2 )m B b^j 
■h (^/x 3 (l-r 2 )m B b 3 ^ + 0(b 3 - h)K (y/x 3 (l - r 2 )m B 6 3 ) 
■h (^3(1-^561)! St(x 3 ), (Al) 



h f (x 1 ,x 2 ,x 3 ,b 1 ,b 2 ) = |6>(6 2 - b 1 )l (m B y/x^x 3 ~{l - r 2 )6i)K (m B a/xiX 3 (1 - r 2 )6 2 ) 

1 / K (m B F (1) 6 2 ), for F 2 > 0\ 
+ (6 ^ 62) }-(fH« (mB ^6 2 ),fo r % <oJ- < A2 > 

where Jo is the Bessel function, K and Jo are the modified Bessel functions with 
K (—ix) = —(tt/2)Y (x) + i(ir/2)J (x), and F?^ is defined by 



(xi - x 2 ){x 3 + (x 2 - x 3 )r 2 ) + rj! 



(A3) 



The threshold resummation form factor S t (xi) can be found in Ref. 1301. 
The Sudakov factors used in the text are defined as 



S ab (t) 



s cd (t) 



s (xim B / V2, 61) + s (x 3 m B / y/2, 63 J + s ((1 - x 3 )m B / y/2, h 
1 

~7i 



In Ht/A \ I In ln(t/A) 



-ln(6 x A) -ln(6 3 A) 
s (xxUIb/ V2, 61) + s (x 2 m B / \/2, 6 2 ) + s Ml - x 2 )m B j a/2, 6 2 

+s (x 3 m B / V2, 61) + s ((1 - x 3 )m B j a/2, 6 
1 



(A4) 



2 In ln(t/A) +ln ln(t/A) 
-111(6! A) + -ln(6 2 A)J 



(A5) 



where the function s(q, 6) are defined in the Appendix A of Ref. [27(]. The scale tj's in the 
above equations are chosen as 



max 



X3 


[1 


— r 2 )m B 


Xi 


[1 


— r 2 )m B 



t 



max 



W Xix 3 (l - r 2 )m B , a/ (x 1 - x 2 )x 3 (l - r 2 ) + r 2 m B , I/61, l/6 2 ), (A6) 



where r = mM 2 / m B{M 2 = J/\I/ ,i] c ),r c = m c /m B . The scale Vs are chosen as the maxi- 
mum energy scale appearing in each diagram in order to kill the large logarithmic radiative 
corrections. 
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APPENDIX B: INPUT PARAMETERS AND WAVE FUNCTIONS 



The masses, decay constants, QCD scale and B meson lifetime are the following 



f J/9 = 0.405GeV, f K = 0.16GeV, f Vc = 0.42GeV, m w = 80.41GeV, 



m r . 



2.98GeV, m B = 5.2794GeV, mj^ = 3.097GeV, 
1.643ps, r B o = 1.53ps. 



For the CKM matrix elements, here we adopt the Wolfenstein parametrization for the 
CKM matrix, and take A = 0.2257, A = 0.814, p = 0.135 and fj = 0.349 Q. 

As for B meson wave function, we make use of the same parameterizations as in 
Ref. [27J. We adopt the model 



!>s(x, 6) = N B x 2 (l — a;) 2 exp 



m B x 



: {u) b bf 



:b2) 



where Uj, is a free parameter and we take cut, = 0.40 ± 0.04 GeV in numerical calculations, 
and N B = 91.745 is the normalization factor for u>b = 0.40. 

For the vector J/ty meson, we take the wave function as follows, 



x 



mj/vt L<j>j/*(x) + i lP<I)j/*(x) \ ■ 



(B3) 



Here, <ft L denote for the twist-2 DA's, and 0* for the twist-3 ones, both of them have 
experimental and theoretical basis (3l|. x represents the momentum fraction of the charm 
quark inside the charmonium. The J/^f meson asymptotic distribution amplitudes read 



as 



a 



X 



0/)8-^Lx(l - x) 



10.94 



2,/2N c ' 



x(l — x) 



2y/2N c 



;i-2x) : 



1 - 2.8x(l - x) 
x(l — x) 



0.7 



2.8x(l - x) 



0.7 



(B4) 



It is easy to see that both the twist-2 and twist-3 DA's vanish at the end points due to 
the factor [x(l — x)] ' 7 . 

For pseudoscalar meson rj c , the wave function is the form of 



(B5) 



The twist-2 and twist-3 asymptotic distribution amplitudes, <p v and s , can be written 



as 



31 



KM 



Vc 



X 



9.58 
1.97 



>lc 



frjc 

2y/2W c 



:X(1 — X 



x(l — x) 



1 - 2.8x(l - x) 

X{1 — X) 

2.8x(l-x) 



0.7 



(B6) 
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The twist-2 kaon distribution amplitude 0^-, and the twist-3 ones and 0^ have 



been parameterized as [32 

fx 



<f>%(x) 



2y/2N c 



6x(l - x) 1 + a? Cl ,2 (t) + a*C*'\t) 
+afC 4 3/2 (t) 
1 + ( 30% - Cl /2 {t) 

3<|^3 + ^(l + 64)j>C< 4 1/2 (t) 



,Kr< 3 / 2 / 



4> T K {x) 



K 



2^2N C 



(1 -2x) 



1 7 3 

1 + 6 ( 5% - -r] 3 uj 3 - — p\ - -p\a 2 



x(l - 10x + lOa; 2 



with t — 2x — 1, the Gegenbauer moments af^ = 0.17, af = 0.115, and af = 
parameters r] 3 = 0.015, u 3 = —3.0, the mass ratio px = (?nu(d) + m s )/mK 
and the Gegenbauer polynomials C^{t), 



Cl'\t) = 3t, C^{t) 



1/2/ 



i(3t 2 -l), Cf(t) = ^(5t 2 -l), 
15 



cl /2 (t) 



(3 - 30t 2 + 35t 4 ) , Cl ,2 (t) = ^ (1 - 14t 2 + 21t 4 ) 



(B7) 



:b8) 



(B9) 

-0.015, the 

: m K /m 0K 



(B10) 
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